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SUPERSONIC FLOW AROUND CIRCULAR CONES AT ANGLES OF ATTACK 1 


By Antonio Feesi 


SUMMARY 

The flow around cones without axial symmetry and moving at 
supersonic velocity is analyzed. Singular points are shown 
to exist in the flow around the cone if no axial symmetry exists. 
The results of the analysis are applied to the determination of 
How around circular cones at small angles of attack. The 
concept of a vortical layer around the cone at small angles of 
attack is introduced , and the correct values of the first-order 
terms of the velocity components are determined. 

The method used is applied to cones at finite angles of attack, 
and it is shown that good agreement with experimental results 
can he obtained from the first-order theory if the complete 
equation for the pressure distribution is used. 

INTRODUCTION 

The flow around a cone having a circular cross section and 
moving at supersonic speed has been determined by means 
of the assumption of small disturbances or by means of 
more rigorous methods that consider the existence of the 
shocks. The latter methods can be applied for any Mach 
number larger than unity and have been developed by several 
authors, at first by assuming all the flow as potential flow 
(references 1 and 2) and later by also considering the variation 
of entropy due to the change in angle of attack (reference 3). 
By means of the development given in reference 3, values of 
flow properties around circular cones at an angle of attack 
have been tabulated in reference 4. The method has been 
extended in reference 5 to larger angles of attack. 

In the method given in references 3, 4, and 5, the flow 
properties were considered continuous and were developed 
in Fourier series in terms of the angle of attack; however, 
the existence of a singular point at the surface of the cone was 
neglected. The derivatives of the flow properties were 
obtained by differentiating the Fourier series term by term, 
and the terms of the series that represent the derivatives 
were assumed to be of the same order as the corresponding 
terms of the integral quantities. For this reason, an errone- 
ous distribution of the entropy at the surface of the cone was 
obtained. 

In this report-, the flow around the cone in the general case 
is discussed, the existence of singular points in the flow is 
proved, and a different procedure for determining the flow 
around cones at small, but finite, angles of attack is devel- 
oped. This procedure shows the way in which the values 
tabulated in reference 4 can be used if a simple correction 
is introduced. The values obtained in this way are compared 
with experimental results at several values of angle of attack. 


SYMBOLS 
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Ti 
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a 


polar coordinates (see fig. 1) 
polar velocity component in radial direction 
(along r), referred to limiting velocity (see 
%• 1) 

polar velocity component normal to v T in meridian 
plane 0= Constant, referred to limiting velocity 

(fig. 1) 

polar velocity component normal to meridian 
plane 0=Constant, referred to limiting velocity 

(fig- 1) 

time 

pressure 

density 

entropy 

ratio of specific heats (cflc,) 

specific heat at constant pressure 
specific heat at constant volume 
speed of sound 

projection of streamline on sphere r=Constant 
with center at center of the polar coordinate 
system 
local velocity 

undisturbed velocity, referred to limiting velocity 
limiting velocity (velocity for expansion in the 
vacuum) 

semiapex angle of conical shock 
inclination of axis of conical shock with respect 
to free-stream direction 

inclination of axis of conical shock with respect to 
axis of body 
angle of attack 


Subscripts: 

1 stream conditions 

a zero-order terms of Fourier series (part independ- 

ent from angle of attack) 

b first-order terms of Fourier series (part propor- 

tional to angle of attack) 

c higher-order terms of Fourier series or quan- 

tities at surface of cone 

e quantities at external surface of vortical layer 

s quantities for polar coordinate system having 

axis coincident with axis of conical shock 


'Supersedes NACA TN 2238, "Sopersonlo Flow around Circular 0 ernes at Angles of Attack” by Antonio rent, 1950. 
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A prime is used to designate the terms of zero order in the 
power series in Ad for the quantities in the neighborhood of 
the meridian plane tf=ir;-two primes - are used to designate 
the factor of the term containing Ad 2 / 2 in the same power 
series. 

THE FLOW FIELD FOR CONICAL FLOW WITHOUT 
AXIAL SYMMETRY 

In order to analyze the flow field for conical flow without 
axial symmetry at supersonic speeds, assume a polar co- 
ordinate system (,r,\p,6). Call v r the velocity component in 
the radial direction, v u the velocity component in the direc- 
tion normal to r in the meridian plane <9= Constant, and w 
the component normal to the meridian plane (fig. 1); that is, 


t>r = 


»«=' 


dr 

dt 

rd\p 

dt 


If the flow is conical, 
dr 

|2_ 0 

dr 


rdd. , 
W ~~df 8m ^ 


bv n „ dio 
-p=0 =0 

dr dr 

f-0 

dr dr 


For these conditions Euler’s equations become 

dv r , w d v r , , „ 

Vn ^rr + -r—i — v n 2 —w I =0 
b\p smf be 


bv„ , w dy„ . 1 dp . . . 

\— V r Vn— ufcoi $= 0 


Vu' 


V * di p 1 sin \p be ^ p b\p 
d w , w bw 


1 |f+tM0+y n w cot^O 


b\p n sin^ d0^p sin^ be 
and the continuity equation becomes 

2 p» r sin^+«nsin^^+psin^^+ VnP cos^+ 

dp , bw . 

W b8 + P b8 = ° 


(la) 

(lb) 

(lc) 


( 2 ) 


Because the energy in the flow is constant, the following 
relations must apply: 

_jY_/lbp p bp\_ ( dy r 

7-iV p 30 P^d) \ Vr be +Vn be + be) ^ 

7 /I dp p bp\ / de,, d»» . bw\ . , . 

7-1 \ P bf p 2 b\p) b* + * bf + bf) (3b ^ 



Combining equations (1), (2), and (3) results in 




+ 


bw 


sin \pbd 


(l-- 

\ WVnf bVn - dw\ 

V a\ 

/ a 2 \sin p be~*~ bp / 


The entropy at any point of the flow can be expressed as 
S R ■ log, -2. f -f Constant 

7—1 b Pi\p/ 

where p and p are the local quantities and pi and p t the 
stream quantities. Therefore, 


7 — 1 bS _ 1 dp 1 dp' 

R be~p be 7 p be 

7—1 bS_ 1 dp 1_ 

R bp p bip 7 p bp J 


(5) 


Combining equations (5) with equations (1), (2), and (3) 
results in the following expressions : 

a 1 bS . . bw bv r bo n . 
y, R sm * bp~ Vr be-** be +VrW sin * +p -* w coa * 


a 2 bS 


bv r bw 


w da* 


yRbp^ Vr bip W bp^ sin p bd 


+*V v a —w* cot p 


Equations (6) combined with equation (la) give 

. , bS bS 
VnSin + —=- w — 


( 6 ) 

( 7 ) 


Equation (7) is general for any conical flow’ and defines the 
lines of constant entropy, which correspond to the stream- 
lines. In fact, if A is the streamline projection on tho sphere 
given by r= Constant, 

dS_^S_dp,^Sde_ a 

dL~bp dL + bd 
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and, from equation (7), 


/ d 6 \ w_ 

WA — ®« sii 


At the surface of any conical body the component of 
velocity in the direction normal to the surface is zero and 
the stream moves tangentially to the body; therefore, the 
entropy at the surface of the body must be constant or must 
change in a discontinuous way (in which case equations (7) 
and ( 8 ) are not valid). 

THE PROPERTIES OF CONICAL FLOW WITHOUT AYTAT, 

SYMMETRY 

In order to analyze the properties of conical flow without 
axial symmetry, consider first a polar coordinate system 
having its axis coincident with the axis of a circular cone 
at a small angle of attack (fig. 2 ) and assume that the direc- 
tion of the ' undisturbed velocity Vi is in the plane 6 = 0 , 
0 =x. In this case, the plane 6=0, 8—t is a plane of sym- 
metry of the flow and, in this plane, 


b6 ° 




^==0 

2>6 


be 


Therefore, equation ( 7 ) shows that in the plane 6=0, 0 =ir 
the entropy is constant. At the surface of the cone (V'=^e), 
the normal component i\ is zero and wt* 0 ; therefore, equation 
(7) shows that the entropy remains constant also along the 
surface of the cone (^=^ c ). Only at points A and B, 
(defined by 0 = 0 , 0 =x, and and w— > 0 ; therefore, 


equation (7) is indeterminate. Because the body is at an 
angle of attack, the axis of conical shock does not remain 
coincident with the direction of the velocity Vi, and the 
entropy in the plane 0=0 must be different from the entropy 
in the plane 0 =x. The entropy at the cone surface therefore 
must be different from the entropy at the plane 0 = 0 , from 
the entropy at the plane 0 =v, or from both. In this case, 
where w=0 and i' a = 0 , a discontinuity of entropy must exist 
either at A or B or both points. 

In order to find a relation between the value of the entropy 
at the surface of the cone and the values of the entropy in 
the meridian plane 0 = 0 , 0 =v, the following considerations 

can be used: In the meridian plane( 0 =O, 0 = r), w, and^y 

are zero because the plane is a plane of symmetry of the flow 


field and, from equation (7), ^-=0. 


Therefore, in the plane 


of symmetry in the zone outside of the singular points A 
and B, 


A /Ag\ 2 da 

b\lf\bd 2 ) sin i (/ b(P 


For the case considered, the velocity component o n is nega- 
tive at the shock or at the Mach cone and remains negative 
throughout the field, until it becomes zero at the surface of 

b 2 S 

the cone; therefore, the value of tends to increase in 

absolute value as $ decreases from the value corresponding 
to the value at the surface of the shock to the value at the 

surface of the cone when ~ is negative and tends to decrease 

, bw . . . 

when -Q 0 is positive. 



Figure 2. — The ringnlar paints st the surface of the cone! 
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Because the entropy remains constant along each stream- 

ers' 

line, the decrease of the absolute value of as ip decreases 


corresponds to a departure of projection of the streamlines 
on the sphere r= Const ant from the plane of symmetry; but> 


is of sign opposite to the projection of the 

streamlines on the sphere r= Constant tends to converge 
toward the plane of symmetry as \p decreases from the value 
at the shock to the value at the surface of the cone. 

Now, with the convention used in figure 2, the component 
w is negative throughout the field and is zero at 0=0 and 

0= jt. Therefore, ~ is negative in the zone 0=0 but is 


positive in the zone 0=ir, and the streamline projection tends 
to converge toward the zone of point A and diverge from 
the zone of point B. Because of the departure of the stream- 
lines from the plane 0=ir, the entropy in zone B remains 
constant and, therefore, the entropy at the surface of the 
cone is equal to the entropy at the meridian plane 0=r; at 
point A a discontinuity of entropy exists from the value 
corresponding to the plane 0=0 to the larger value existing 
in the plane 0=ir. All the projections of the streamlines 
converge at point A where the entropy is not single-valued. 
Because v n approaches zero near the cone, equation (8) 
shows that all the streamline projections tend to become 
parallel to the line ^= Constant in the zone near the cone 
and converge at A. The value of v n in the meridian plane 
0=0 can change sign and can be positive in the neighbor- 
hood of the point A (case of large angles of attack). In this 
case the right-hand side of equation (9) changes sign and 
the singular point moves away from A in the meridian plane 
0=0 and occurs at the other point of the meridian plane 
where v n is also zero. (At the shock or at the Mach wave 
v n is negative; therefore, another singular point where v„=0 
must exist.) 

It is interesting to observe that singular points must exist 
in any supersonic conical flow without axial symmetry. 
Considerations similar to those used for the case of circular 
cones at an angle of attack can be extended to other cases, 
and it can be shown that the streamlines that are tangent 
to the body start from points of the shock and meet the 
body at points where the component, of velocity perpen- 
dicular to the radius and tangent to the body vanishes and 


has a positive derivative in the streamline direction ^equiva- 
lent to the condition of positive Convergency of 

streamlines occurs and, therefore, the points are singular at 
the points where this component vanishes and has a negative 
derivative in the direction tangent to the body, while the 
component normal to the body also vanishes and has . a 
negative derivative in the direction normal to the body 

^equivalent to v n =0 and negative^- 

For example, the conical body of figure 3 has two planes 
of symmetry, AA' and BB' when w and v n are zero, but at 

the points BB', is positive, while at AA', it is negative. 


Therefore, because ^ is negative at AA' and BB', the 

points AA' are singular points and the entropy at the surface 
of the body is determined by the shock strength at the 
points CC'. 

DETERMINATION OF THE FLOW AROUND CIRCULAR CONES 
AT SMALL ANGLES OF ATTACK 

In order to determine the flow around circular cones at 
small angles of attack, consider a polar coordinate system, 
the axis of which is coincident with the body axis. At the 
surface of the body the velocity component v K is zero and in 
the neighborhood of the body is very small; therefore, the 
terms v^/a? can be neglected with respect to unity. 

If the angle of attack is small, the component w is also 
small and the terms vf/a 1 can also be neglected. On the 
basis of this approximation, in the neighborhood of the body 
equation (4) can be expressed as 

2„,+„. cot dO) 


This equation permits a particular solution of this type 
chosen from physical considerations 


V T = V r + av fi COS C0S m 0 

cos 8+'Y',o._ cos mO - 
w—aw b sin 0+S w c cos md 


(ID 


where v u b t , v n ^ be , and w„_ c are functions only of ip, are 
constant for constant values of ip, and must be chosen in a 
form that satisfies the boundary conditions. 

Consider now a conical shock having circular cross section 
and semiapex angle a (fig. 4). Consider the cone inclined at 
an angle S with respect to the undisturbed velocity with a 
polar coordinate system, the axis of which is coincident with 
the axis of the conical shock. If d*,, and w, are the 
velocity components referred to the limiting velocity I i in 
the new coordinate system (r„ ip„ 0,), from the shock-wave 
relations the following equations result: 


v, t = Vi cos S cos o- + Vi sin 5 sin a cos 6, 


w,= — Vi sin 5 sin 0, 


—v a 


y—l 1 — g,/— W,* 

y+1 Vi cos 5 sin <x— W sin 8 cos a cos 0, 


If 5 is assumed to be small and terms of the order of 5* are 
neglected, these equations become 


v, = Yi cos <r+ 5 Vi sin <r cos 0, 

w, =— Fi* sin 6, 

Y— 1 /l — Fi s cos 2 ff\ 
r,t ' - 7+1 V Fisin o- } 

1=1 ! cos , (- 2 F, + 1 -y-‘. c « l£) cos ). 

7+1 V Fi sm* a / 


( 12 ) 
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F:pcbe 3— The atngiilar points at the surface of* cone in supersonic (low without axial symmetry. 



Florae 4. — The position of the conical shock with respect to the conical body. 


The flow behind a circular conical shock inclined at a small 
angle 5 with respect to the undisturbed stream can therefore 
be expressed in the form 


COS e, 

Vn,= V*a,+ Sv *l, s <iOSe * " 

w, = Svjt, s sin 8, „ 


(13) 


where all the terms containing v n , v„ b , and w 6 are small, 
that is, of the order of 5. 

If the axis of the conical coordinates is rotated at an angle 
V of the same order as the angle 5 and terms of the order of S 2 


are neglected, the velocity components v T , v„, and w referred 
to the new axis (fig. 4) become 



where 8, the coordinate referred to the new axis, is given by 
cos 0=cos 8,— i] cot sin 2 8, 

Equations (14) show that the flow behind a circular 
conieal.shock inclined at a small angle 5 with respect to the 
direction of the undisturbed stream and at an angle y with 
respect to a conical body can also be expressed in the form 
given in equations (11) where, if the angles y and 8 are 
small and the terms of the order of y 2 and 5 ! are neglected, 
only the terms having the subscripts a and b must be con- 
sidered. Therefore, equations (11), when the terms with 
the subscript c are neglected, are valid for small angles of 
attack ( 01 = 5 — 17 ), and a conical circular shock is consistent 
with the solution chosen for the flow around the body. At 
the surface of the cone the assumption that to* is small 
corresponds to the assumption that only the first-order terms 
of angle of attack are considered. The conical shock is 
inclined at an angle 77 with respect to the circular cone. 

This analysis is similar to the analysis of references 2 and 
3. No assumption, however, has been introduced for the 
entropy distribution ; only the velocity components have been 
considered to be in the form of equations (11), and no 
limit ing assumption has been introduced for the derivatives. 
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In reference .3, in addition to the equations 

Vr=Vr a +CtV ri COS 9 (15a) 

v n =v H +ocv ni cos 0 (15b) 

w=aw s sm.5 (15c) 

the expressions 

p=p a +ap b cos 9 (16a) 

p = p a + aps COS 0 . (16b) 

S=S a -\- aS b cos. 0 (16c) 


have been used, and the derivatives of entropy, pressure, 
and density have been obtained from differentiation of 
equations (16). In this way a solution has been found which 
gives values of entropy that are variable along the cone 
surface and are constant in each meridian plane, while the 
entropy actually remains constant along the cone surface 
and changes in the meridian plane. An incorrect entropy 
distribution has therefore been obtained at the surface of 
the body. 

In order to analyze in more detail the significance of equa- 
tions (16) and their inconsistency with the approximation 
considered in references 2 and 3, consider equation (7). In 

bS 

the plane of symmetry 0=0 or 0=r, ^=0; therefore, 

S tt +aS b or S a —aS„ of equation (16c) must remain constant. 

Consider now the plane 0=ir and express the entropy S 
in the form 

S=S'^S"^f • (17) 


which satisfies the condition of symmetry. Because of 
equation (7), 

bS' 


bi 


=0 


and, from equations (16) and (17), 

S'=S a -aS b 


and 


aS b =—S n 


From equations (7) and (17), 

'i)S' bS" A6 


( 18 ) 

( 19 ) 

( 20 ) 


. , /bS' bS" A 


However, *^-=0; therefore, 


Ad bS" 
2 b$ 


wS” 

v n sin ^ 


or, since from equation (15c) w= —aw b Ad, 

bS" aWtS” 

bip sin f 


By use of equation (20), in the neighborhood of the plane 

9 = T, 

bS b 2w b S b a i 2aw b S" 

a -5-7-= = 7 = : 7 

o\j/ v n sm 5,sm ^ 


In reference 3, the term 


wS " 


has been considered every- 


th sm \f/ 

where to be of the order of a ! and has been neglected ; hence, 


bS b . , bS a 

-5rr= 0 and ■x-r=0- 


bip 


b\(> 


However, neglecting this term is cor- 


rect only when the ratio Wt/v n is of an order different from 
1/a and, therefore, when |a«|»0. Near the surface of the 

cone, v n approaches zero and, therefore, the term be- 
comes large and cannot be neglected. The extent of the 

bS" 

field around the cone where the term is of the order of 

0\{/ 


a can be easily determined. 

Consider a polar coordinate system having its axis coin- 
cident with the axis of the cone. At the surface of the body 
p tt is zero; therefore, in the neighborhood of the body the 
velocity component v n can be expressed in the form 



or, by use of equation (10), in the form 


( 2 »,+ s1 ^| 55 )(4«k (21) 


Therefore, v n is of the order of a when (A^) r is of the order 
of a. InJhis conical layer of thickness (A yf/) T of the order of 

is also of the order of a or larger. In this layer, which 


bip 


bS . 


can be called the vortical layer, the term also is of the 
order of a because from equations (18) and (19) can be 


shown to be of the same order as 


bS " 
bf 


In order to investigate the effect of this vortical layer on 
the velocity and pressure distribution at the surface of the 
cone, consider equations (2), (4), (6), and (7). If the velocity 
components in the neighborhood of the meridian plane 9=r 
are expressed as 


, A 9* „ 

Vr =V t 2" Vr 

(22a) 

, A 0* „ 

V n = Vn 

(22b) 

s 

II 

s 

-s 

1 

«|& 

(22c) 
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the following expressions can be obtained: 


bv/ 


-=v a 


bS " 2 w"S" 

bp v a ' sin p 

2w" 


\Sin p / 


bw' 


bp v n ' sin p 


bp ““ ®*' Vsin P 

(v/v" + v n 'v n " + w"v/ sin p+v a 


’w" COS p- 


1 

2 B 

fo n 'v n " A 

— V n " COt 1 p — B r " 


w" 

/2w" 2 2v n 'v„" . 2 Bn sin p bw"' S 

2v/w" 2 l 

bp . Vn' 2 

\ a ' 2 A a' 2 ) 

K 2 a' 2 ) 

sin p 

1. a ' 2 a ’ 2 1 a ' 2 bp J 

a' 2 J 


where 


— v/t l/' — Vn'v/' + w " 2 y — 1 


Fr'+V a ' COt p- 


V) 


B-- 


sin p 


1 -- 


and 


bf 


=-v.— 


a' 2 1 bS a 
Vr. yR bp 


(27) 


(28) 


(29) 


Therefore, all the derivatives of the velocity components 
of zero and first order are affected by the entropy variation 

in the meridian plane.. The terms and-^— are of the 

order of a. in a layer of thickness a near the surface of the 

cone. In this layer they change the value of of quantities 

of the order of a and the pressure at the surface of the cone 
of quantities of the order 8 1 or. The only place where the 

terms and 3 - 7 - are of zero order is near the surface of 
bp bp 

the cone in a layer of thickness a 2 , where v n tends to the order 
0 /a* (equation (9)). Because the effect of this vorticose 
layer existing at the surface of the cone on the pressure is 
of the order 0 /a 2 , it can be considered in this approxima- 
tion, which neglects terms of the order of a 2 or higher, 
that across this layer the pressure distribution remains 
constant, but an abrupt variation of entropy occurs; there- 
fore, in this approximation the phenomenon can be repre- 
sented as in references 3 and 4, where the entropy remains 
constant in every meridian plane until a vortical layer of 
infinitesimal thickness is reached at the surface of the cone 
across which a variation of entropy occurs from the value 
5,+ aS,, cos 6 to the value S a —ocS t that exists at the surface 
of the cone. Across the layer a variation of density and 
velocity components occurs and can be easily determined. 

Let F re , w„ S e , p„ p„ and a* be the quantities at the external 
surface of the layer (these are the quantities tabulated in 


(23) 

(24) 

(25) 

(26) 


reference 4) and v Tc , w c> S c , p c , and p e be the quantities at the 
surface of the cone. Because it has been shown that 


S'.' a' 
yR 




then 


P'=P* 


Vr e 2 + W e 2 — V r 2 — W c 2 = 


y—1 


a. 


S'- Sc 


(30) 

( 31 ) 


where S,—S e is the entropy jump across the layer. 
Now, 

v r =v Ta +av r 6 cos 8 
w=aw t sinff 
and in the plane 8= tt 

»r.=v rr _ 


w e =w c —b 

S e =S e 

(»r e -«lV 4 ) e =(«r e -aPr 6 ) c =V/1 

(S'-aS t ) =(£„-*£>) =S' J 


or 


(32) 


Therefore, if terms of the order of a 2 are neglected, equation 
(31) becomes 




a 2 ^ 

y—1 c P 


( 33 ) 


or 


a- 


S, 


However, from equation (la), at the surface of the cone, 

bv r 


or 


w= 


w t = 


sin p b6 


sin p e 


(35) 


therefore, the values of tv and w b can be determined from 
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equation (34) and from the tables of reference 4 where 


S>=C,(%- 7^) 

(36) 

\Pa Pa/ 

a t 2 = y ~-{ 1-r/ 2 ) 

(37) 


The values v r<l , v, b , and w b having been determined, the values 
of v, and w can be obtained at any part of the cone from 
equations (15). 

THE NUMERICAL DETERMINATION OF THE FLOW FIELD 
AROUND CONES AT SMALL ANGLES OF ATTACK 

The method presented permits the determination of the 
pressure distribution around the cones with the assumption 
that the terms of the order of can_.be neglected. The 
pressure at any point can be obtained from the equation 


£ 

Pi 


' 1 — F 2 Y=i - 

J=V?) e 


1 

7“1 


AS 


(38) 


where F is the local velocity corresponding to the pressure p 
and AS is the increase of entropy with respect to the stream 
conditions where pi and Fi exist. From equations (15), 
F 2 is found to be 


V‘=v r 2 +cc l v, 2 cos 1 8+2aVr a v ri cos 0+ ahoj 2 sin 2 0 <39) 

However, in equation (39) the terms of the order of a 2 
having the form 2v ra v, e cos md have been neglected (see 
equations (11)). If all the terms of the order of a 2 are 
neglected, V 2 becomes 

V i =V, a * + 2aVr a Vr t COS8 (40) 

Equations (39) and (40) are different in terms of the order 
of a 2 ; however, for finite angles of attack good agreement is 
obtained only if equation (39) in which some of the a 2 
terms are retained is used in equation (38). .. 

The reason for the better approximation given by equation 
(39) can be understood if the magnitude of the terms con- 
taining c? in the expression of V 2 is considered. Along the 
surface of the cone, aw t is given by 


aw* 


sin^ 


(41) 


only term retained in equation (40). Because the term 
c?w s is correct and significant, it can still be retained also 
when other terms in a? are neglected and, therefore, equa- 
tion (39) is the expression that must be used for finite angle 


of attack. ^For example, fora 10° cone (^=10°) sinV c =-^ 


and ;%: == 1 for <*=1.75°, which can be considered a small 
angle.^ 

In reference 3, equation (16a) has been used in the deriva- 
tion of the method; however, the use of this equation is not 
necessary, as is briefly shown in the following paragraph. 

Consider a conical circular shock having the semiapex 
angle <r and inclined at an angle 6 with respect to the undis- 
turbed stream. In the neighborhood of the plane 0=*- the 
velocity components can be expressed as in equations (22) 
and at the shock are given by (from equations (12)) 




v/—Vi cos ar—\\ S sin a 
v,"= — fiFi sin cr 
w"=T \S 
y— 1 1 — Fi* cos 2 er 


7 + 1 Fi sin a 
7 — 


7 + 1 

//■ 7—1 
* 7 + 1 


(42) 

(43) 

(44) 


8 cos v ^ 

-2FH 

1 — TV COS 2 cr\ 
Fi sin 2 <r ) 

(45) 

a COS <r ^ 

-2 TV 

, 1 — Fi* cos 2 o-\ 
Fi sin* rr ) 

(46) 


and S' and S" are given by 

S'=c, log, [ 1 ■ +( L ~ M, 2 + 1) [TV sin 2 (<r+ a)— rV 2 ! } + 
{c v -c,) log, T 7 l8in (‘ ff+4 ) - (47) 

TV cos tr sin a—Vn'— t 
S"=28 


l+Xl Ml t 


Fi 2 sin 2 o-— o a /2 


(48) 


For a finite value of a and a small value. of &, is of 

the same order as av Ti because sinV, is also small; therefore, 
c+Oj 2 can have an effect on the velocity and pressure dis- 
tribution of the same order as the term aP, a iv s which is the 


All the velocity components are referred to the limiting 
velocity, and F t is the undisturbed velocity also referred to 
the limiting velocity. 

In the meridian plane 0 =t the entropy S' is constant; 
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therefore when tr=0, from equation (la), 



If equation (4) is applied to the meridian plane 9=t and 
equation (21) is used, the radius of the hodograph diagram 
at any value of ^ smaller than a is given by (reference 6) 


values of reference 4 for the conditions outside the vortical 
layer, and the pressure at the surface of the cone has been 
determined in the following way: 

From tables of reference 4 the value of 5/a has been deter- 
mined (5/e of tables). The position # in the plane B—v of 
the conical body in the shock coordinate system is 

\{f=\pc+ct—S (54) 


C B)f= — 


Br'+Vn COt fahjZ 


W 


sm ^ 


where 


7 — 1 


l- : 


a— 1 _v r —L\ 


(50) 


(51) 


Therefore, at the point f—A $ in the meridian plane 8 = v, 


The value of v Ta at fa-\-a—8 has been determined from the 
tables of reference i(v r at ^ is given by u s —2u 3 ^ ^ 


*> 


The value of r" has been obtained from the tables in shock 
coordinates 

7T- [**.+("0* 3 <55) 


(®r')*- A* =—(»!.')* sin A f + (v r '—R)+ COS A lH-(B)* (52) 

' and 


(.Vn%-M=(Vn')t cos A^+(o/— R)+ sin A<]/ (53) 


From equations (50) to (53), the velocity components c/ 
and v n ' can be determined if the component vf is known in 


the meridian plane 6— -nr. 


Since 


bw" 


is given by equation 


(25), the value of w" can be determined for any value of fa 
But v Ta and v„ tt are the quantities obtained for zero angle of 
attack at the same coordinate ^ of a coordinate system in 
axis with the conical shock, and 


r r n =v/-i\ 


r •—.« ' — r 
t n — * n 


therefore, the entire flow field can be determined until i\' 
becomes 0. The value of ^ for which v a '=Q corresponds to 
— 5 and, therefore, gives the value of a. 

Considerations similar to those used for cones can be used 
for the characteristics method presented in reference 6. In 
this case the pressure at the surface of the body can be 
obtained from the complete equations, and the vortical 
layer must be considered in order to obtain the correct distri- 
bution of entropy. The application is the same because the 
entropy does not change along each streamline. 

COMPAEISON WITH EXPERIMENTAL RESULTS 

In order to have an indication of the accuracy that can be 
expected from the first-order theory, theoretical results have 
been compared with some experimental results available. 
The theoretical results have been obtained by using the 


Then (®r a ) c , (ur 5 ) e , and (w? 6 ) e in the body coordinate sys- 
tem have been obtained by means of equations (32), (34), 
and (35). 

From equation (36) S„ has been obtained (it has a negative 
value), and from tables of reference 7 S a can be determined — 
for example, from the value of the angle of the shock obtained 
from reference 7 — and 


AS=S-aS t (56) 

The pressure has been obtained from equation (38). 

In figure 5 a comparison is presented for a cone of fa= 7.5° 
at M— 1 .6 and four angles of attack. The experimental data 
are obtained from tests performed in the Langley 4- by 4-foot 
supersonic tunnel. For comparison the values given by 
reference 4 and by linear theory are also shown. In figure 6 
a comparison is presented for a cone of ^=10° at M=6.86 
and two angles of attack. The experimental data have been 
obtained from tests performed at the Langley 11-inch hyper- 
sonic tunnel. The agreement in both cases is good, even at 
angles of attack where it would be expected that higher-order 
terms would be important. 

CONCLUDING REMARKS 

The flow around cones -without axial symmetry at super- 
sonic velocity has been analyzed. Singular points which 
complicate the analysis of the flow field were shown to exist 
in the flow. The results of the analysis were applied to the 
determination pf the flow around circular cones at an angle 
of attack. The concept of a vortical layer around the cone 
at small angles of attack has been introduced, and the correct 
values of the first-order terms of the velocity components 
were determined. 



Pressure coefficient 
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FifSCBK 6. — Comparison ot experimental pressure distribution over the soriece of a cone with the theoretical pressure distribution. ^-.—10°; Af— BJifi. (Experimental values for two different 

values of r along the cone surface; data obtained from the Langley U-Inch hypersonic tunnel.) 


The method determined was applied to cones at finite 
angle of attack, and it is shown that good agreement with 
experimental results can be obtained from the first-order 
theory if the complete equation for the pressure distribution 
is used. The analysis can be extended to the application of 
the characteristics method around bodies of revolution at 
small angles of attack. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ya., September 26, 1960. 
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